Out-of-plane thermopower of strongly correlated layered systems: an application to 
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We calculate the out-of-plane thermopower in a quasi-two dimensional system, and argue that this 
quantity is an effective probe of the asymmetry of the single-particle spectral function. We find that 
the temperature and doping dependence of the out-of-plane thermopower in Bi2(Sr,La)2CaCu20g+i 
single crystals is broadly consistent with the behavior of the spectral function determined from 
ARPES and tunneling experiments. We also investigate the relationship between out-of-plane ther- 
mopower and entropy in a quasi-two dimensional material. We present experimental evidence that 
at moderate temperatures, there is a qualitative correspondence between the out-of-plane ther- 
mopower in Bi2 (Sr,La)2CaCu208-t-«, and the entropy obtained from specific heat measurements. 
Finally, we argue that the derivative of the entropy with respect to particle number may be the 
more appropriate quantity to compare with the thermopower, rather than the entropy per particle. 

PACS numbers: 74.72.-h, 72.15.Jf, 65.40.gd 



I. INTRODUCTION 

As is well known, high-temperature superconductors 
(HTSC) share the Cu02 plane as a common structural 
unit, and exhibit highly two-dimensional electronic prop- 
erties This extreme anisotropy is demonstrated by 
the ratio of the in-plane (a-6-axes) and out-of-plane (c- 
axis) resistivities, pd Pah, which can be of the order of 
10^-10^. The unusual superconducting and normal-state 
properties of these materials are thought to originate 
from strong correlation effects in the a-h plane. Hence, 
the in-planc charge transport in these materials has been 
studied much more extensively than the out-otjDlane. 
However, it has been noted by many authors 0, IM Hi 
that the out-of-plane transport can be used to provide 
useful insights about the in-plane physics. Under cer- 
tain assumptions, out-of-plane transport properties de- 
pend only on the in-plane single-particle spectral func- 
tion, and the out-of-plane hopping integral. To date, 
attention has primarily focused on the out-of-plane resis- 
tivity and magnetoresistance, with much less attention 
given to thermal and thermoelectric properties. 

The significance of thermopower has, however been 
recognized in recent years. First, it is a quantity inde- 
pendent of sample dimension, and, like the Hall effect, is 
often insensitive to the grain boundary and/or disorder. 
Second, most of the transition-metal oxides rarely show 
complicated non-equilibrium phenomena such as phonon 
drag, because of poor mobility. Recent work by Kontani 
has suggested that the in-plane thermopower of HTSC 
can be rather simply understood within spin-fluctuation 
theory. 

The electronic contribution to the thermopower of a 
material is fundamentally related to the difference in the 
response of electrons and holes to an applied tempera- 
ture gradient Q ■ The thermopower therefore probes the 
particle-hole asymmetry of the system. The c-axis ther- 



mopower is particularly useful in this regard, since it can 
be directly related to the particle-hole asymmetry of the 
single-particle spectral function. 

In the first part of this paper we will give an explicit 
calculation of the c-axis thermopower within the tun- 
neling Hamiltonian formalism, and show that the re- 
sult can be written as a Mott formula. Using recent 
experimental measurements of the c-axis thermopower 
of Bi2Sr2CaCu208+i (Bi2212) single crystals, we then 
investigate what can be inferred about the asymme- 
try of the single particle spectral function around the 
chemical potential. Together with the c-axis resistivity, 
which primarily reflects the symmetric part of the spec- 
tral function, these two quantities suggest that the spec- 
tral function becomes weakly temperature dependent at 
high temperatures. At lower temperatures, the c-axis 
thermopower may reflect both the physics of the pseu- 
dogap, and superconducting fluctuations. We show that 
the properties of the spectral function inferred from c- 
axis transport are in broad agreement with information 
provided from ARPES and tunneling spectroscopy. 

We note that the opposite procedure, starting with 
the electronic structure, is often used to investigate the 
thermopower. Within a Boltzmann transport picture, 
and under the assumption that the relaxation time is an 
energy independent constant, the thermopower is deter- 
mined entirely by the band structure of the material. 
This approach has been used, for example, to calcu- 
late the in-plane thermopower of the layered cobalt ox- 
ide Naa;Co204. The necessary quantities are either ob- 
tained from electronic structure calculations [8], or by 
using photoemission spectroscopy to determine the den- 
sity of states, combined with some assumptions about 
the energy dependence of the group velocity 0]. 

The second part of this paper explores the relation- 
ship between the c-axis thermopower and the entropy 
of the system. The basis for this relationship are the 
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Kelvin/Onsager relations [13, Ulii which relate the See- 
beck and Peltier coefficients of a material. This will 
be discussed further in section |Vl here, we will use an 
empirical relationship recently noted by Behnia et al. 
p.2i] as motivation. For non-interacting electrons with 
a quadratic dispersion and an energy-independent re- 
laxation time, it can be shown that the thermopower 
S is just the equilibrium entropy S per electron, S = 
—S/NeC. Behnia et al have considered the quantity 
q = eS/C, where C is the electronic specific heat per 
electron/hole (for non-interacting electrons at low tem- 
peratures, C is also equal to the entropy per charged 
carrier). The quantity q is equal to —1 for free electrons, 
and +1 for free holes. Empirically, it was noted that if q 
is calculated for a range of correlated materials from the 
low temperature {i.e. T — > 0) thermopower and specific 
heat, values of q were still typically close to ±1. Miyake 
and Kohno [Ts'l have shown that this universal value of 
q in the T — > limit can be understood within a Fermi- 
liquid picture. Zlatic et al [13] have also investigated 
this relationship in Ce, Eu and Yb-based heavy Fermion 
compounds via a DMFT approach. 

Another limit in which the thermopower is directly re- 
lated to the equilibrium entropy is in finite bandwidth 
systems at temperatures much larger than the kinetic 
energy scale. In this limit the thermopower is related 
to the derivative of the entropy with respect to particle 
number, rather than to the entropy per particle. The 
resulting expressions are referred to as Heike's formulae 

[a [3. 

Despite the existence of 'exact' results in the particu- 
lar cases indicated above, it is clear that the relationship 
between thermopower and equilibrium entropy cannot be 
truly universal. This fact is obviously demonstrated by 
the existence of materials like Bi2212 where the ther- 
mopower is anisotropic. If the thermopower in one di- 
rection is closely related to the equilibrium entropy, then 
the thermopower in other directions will not be. Nev- 
ertheless, given the general difhculties with interpreting 
the thermopower, any relationships which can be estab- 
lished with more familiar quantities are very useful. It is 
therefore important to identify those special cases where 
relationships with the entropy hold, and to understand 
the reasons for this. 

In this paper we will discuss the relationship between 
thermopower and entropy in general terms, using the 
thermopower (calculated in the relaxation time approx- 
imation) and entropy of electrons in a 2d tight-binding 
band as a simple example. In particular, we argue that 
the derivative of the entropy with respect to particle 
number may be the more fundamental quantity to com- 
pare with the thermopower, something which has not 
been previously recognized. We argue that the c-axis 
thermopower of a quasi-two dimensional material is one 
of the special cases where a close relationship with the 
equilibrium entropy might be expected to hold. We will 
see experimentally that at moderate temperatures there 
is a qualitative correspondence with the entropy per par- 



ticle, and that theoretically a similar qualitative corre- 
spondence can be established with the derivative of the 
entropy with respect to particle number. 

The rest of this paper is organized as follows. In Sec. |TT] 
we analyze experimental measurements of the c-axis ther- 
mopower in Bi2212 single crystals from Rcf. The 
calculation of the c-axis thermopower is given in Sec. lIIII 
In Sec. IIVI we discuss the form of the tunneling matrix 
element in Bi2212, and compare the information about 
the asymmetry of the spectral function provided by the c- 
axis thermopower with that obtained from ARPES and 
tunneling spectroscopy. In Sec. |Vl we discuss in gen- 
eral terms the relationship between thermopower and en- 
tropy, using the non-interacting 2d tight-binding model 
as an example. Finally, in Sec. lVII we compare the exper- 
imentally measured thermopower with the entropy from 
specific heat measurements, and theoretically compare 
the c-axis thermopower with the derivative of the entropy 
with respect to particle number. 



II. C-AXIS THERMOPOWER OF BI2212 

Here we consider the c-axis thermopower data from 
Bi2212 single crystals reported in Ref. [l7| that provides 
motivation for this work. Fig. [ija) shows the tempera- 
ture dependence of the c-axis thermopower (S'c) in Bi2212 
with various oxygen contents. The respective values of 
Tc for p = 0.13, 0.16 and 0.20 are 85K, 89K and 84K. 
The carrier concentration per Cu {jp) is estimated from 
the empirical relation to the room-temperature in-plane 
thermopower Sah [3- Unlike Sah [l7], the sign of Sc 
is positive over the measured temperature range for all 
dopings. The key features of the data are i) an approxi- 
mately linear temperature dependence at moderate tem- 
peratures, becoming sublinear at higher temperatures; ii) 
the slope of the linear part is positive, and has a magni- 
tude which decreases with doping; iii) an upturn in the 
thermopower at lower temperatures, with the tempera- 
ture at which the upturn occurs decreasing with doping. 

Having identified the salient features of the data, we 
now turn to the theoretical calculation of the c-axis ther- 
mopower. 

III. THEORY OF C-AXIS THERMOPOWER 

We consider the c-axis thermopower of a quasi-two di- 
mensional material consisting of well-separated planes of 
atoms in which electrons move relatively easily, but in- 
terplanc motion is suppressed. This has two important 
consequences. First, the interplane hopping integral, tc is 
much smaller than the typical in-plane hopping integral 
tab- Second, interactions between electrons in neighbor- 
ing planes are expected to be weaker than between elec- 
trons in the same plane. Further, we will assume that 
each layer is translationally invariant, and hence the in- 
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FIG. 1: Out-of-plane thermopower of Bi2212 single crystals 
(reported in Ref. for three different doping levels, and 

in the natural unit of fcs/e. The doping level is quantified 
by p, the carrier concentration per Cu. Also shown are fits 
to the thermopower over the appropriate temperature range 
(see main text). The inset shows the fit parameters for each 
doping. 



plane momentum is conserved in the tunneling process. 
This assumption could be violated if, for example, the c- 
axis transport occurred via the resonant tunneling of the 
electron through impurities located between the layers 

The most suitable formalism in which to calculate the 
c-axis thermopower is the tunneling Hamiltonian [20j| . In 
this approach, transport quantities are calculated pertur- 
batively, with the small parameter taken to be tc, rather 
than the magnitude of the applied field. Further, if inter- 
actions between the layers are neglected, then there are 
no vertex corrections to the elementary bubble. We stress 
that in using the tunneling Hamiltonian approach with 
tc treated as a small quantity, we have not assumed that 
the c-axis transport is necessarily incoherent Q • The co- 
herence issue depends on the relative size of the Fermi 
energy ep and the inelastic scattering time r [2l[. Since 
ep will depend on both tab and t^ it is completely pos- 
sible to have coherent transport with epT ^ h, despite 
having t^r <^ h. 

We start with the Hamiltonian of the complete system 
(i.e. all layers) in a real space representation, and then 
transform into momentum space along the in-plane di- 
rections only. The Hamiltonian is then written in terms 
of operators aj^ - (we have suppressed spin indices here) , 
which create an electron in the ith layer with in-plane 
momentum k. We assume translational invariance in the 
in-plane direction, and hence the tunneling matrix ele- 
ment satisfies Tpk = tk<^pk- The form of is arbitrary 
at this stage. We then separate out the part of the ki- 
netic energy which hops particles in the c-axis direction 
i.e. transfers particles between layers. We further assume 



that this kinetic term is the only term in the Hamilto- 
nian which couples one layer to the other. This leads us 
to the following Hamiltonian 



H = '^Hi + <ka^,_,_iak,i + h.c. 



(1) 



i,k 



where Hi refers to the in-plane Hamiltonian, which con- 
tains only operators belonging to the zth layer. 

To convert this into a tunneling problem, we argue that 
to lowest order in fk, the current of electrons along the 
c-axis is limited by the current passing between just two 
planes. We can then simply isolate any two planes, and 
use the tunneling approach to calculate the current to 
lowest order in tk- We therefore work with the Hamilto- 
nian 



H' = Hi + Hi+i + ^ tka^i+iOk.i + h.c. = Hq 



Hj 



(2) 



Since the Hamiltonians Hi and i/i+i contain only oper- 
ators in the zth and (« -I- l)th layers respectively, they 
commute with each other and effectively represent inde- 
pendent systems. 

To calculate the thermopower in the tunneling Hamil- 
tonian formalism, we exploit the Kelvin relation between 
the thermopower and Peltier coefficient, H — TS. Re- 
call that H is the coefficient of proportionality between 
the heat current and the electrical current flowing under 
isothermal conditions. The total heat current is related 
to the energy current hy Q = U + {fi/e)I. We there- 
fore use the tunneling Hamiltonian approach to calculate 
the energy and electrical currents between the layers in 
response to an applied potential difference V. The elec- 
trical and energy currents are given by 



e • • ze 



= ie^(ika[ ,+iak||,j - h.c), 



(3) 



U{t) = \m - m+i)) = '-{[Ht,H, - H,+,]). (4) 

Unlike the electrical current, an explicit form cannot be 
obtained for the energy current, since we have not spec- 
ified the form of the in-plane Hamiltonian. However, 
within the tunneling Hamiltonian approach this does not 
matter, since we can obtain a general expression for the 
energy current without specifying the in-plane Hamilto- 
nian explicitly. 

The electrical and energy currents are calculated using 
the Kubo formula with Ht treated as the perturbation. 
We find, to lowest order in ti^: 



-2eIm[H, 



(5) 
(6) 



U{V) = 2lu,[nret{-eV)] - '^Itl^i^ 

le 

H(zc^„) = ^|ikP^^e(k,ze„)g(k,ie„-|-ic^„), (7) 
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FIG. 2: The diagrams evaluated in the calculation of the ther- 
mopower, where ie„ and iujn are the fermionic and bosonic 
Matsubara frequencies respectively. Circles indicate vertices 
proportional to tk, while squares indicate vertices propor- 
tional to fk(*en + iiOn/2). (a) Tl{iujn). (b) f2(iLj„). 



k ^ ie„ 

X C;(k, ie„)fj(k, ie„ -|-ia;„). (8) 

These correlation functions are shown diagramniatically 
in Fig. [21 The absence of interactions between the layers 
means that there are no vertex corrections to these ele- 
mentary bubbles. The thermopower is obtained via the 
Kelvin relation as Sc = limy^o(l/r)Q(F)//(t/): 



as 

Y^tl f &(2A„(k,.)^(k,.))i^(-^) 

^ _ I k 

' J2tlfde{Al{Ke)+Al{k,e))(-^) 

k 

(11) 

In particular, if the spectral function is particle-hole sym- 
metric for all k {i.e. Aa(k, e) = 0), the numerator of 
(|lip vanishes, and the c-axis thermopower is zero. Note 
that the c-axis conductivity, which is (up to a constant of 
proportionality) the denominator of Eq. (fTTj) . also probes 
both the symmetric and antisymmetric parts of the spec- 
tral function. However, because A(k, e) > 0, and hence 
|^s(k, e)| > |Aa(k, e)|, it follows that the conductivity is 
dominated by the symmetric part of the spectral func- 
tion. 

In the next section, we will apply these ideas to Bi2212. 
We will investigate what information about the symme- 
try properties of the spectral function can be inferred 
from the c-axis thermopower, and whether this is consis- 
tent with ARPES and tunneling spectroscopy. 



_ . i:4/^.^-(M)i^(-g) 

^ j:^lJdeA^iKe){-'^) 

k 

where ^(k, e) = — (l/7r)Im[G'^(k, e)] is the single-particle 
spectral function, /(e) = (1 -I- e^'^'^~'^^)~^ is the Fermi 
function, and fi is the chemical potential. This re- 
sult coincides with that found by Nagaosa Q. If 
we define an energy-dependent conductivity crc(e) = 
2e'^J2]^tlA'^{k,e), then the c-axis thermopower may be 
written as 



J deCTc(e)( ^ ^ 
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de) 
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(10) 



which is known as a Mott formula. Note that this re- 
sult has been derived for an arbitrary in-plane spectral 
function, and no assumptions have been made about the 
nature of the in-plane physics. We have shown explic- 
itly that within the tunneling Hamiltonian approach, the 
thermopower may always be expressed as a Mott formula 

Eq. ^ shows that the c-axis thermopower probes the 
asymmetry of the spectral function around the Fermi en- 
ergy. To demonstrate this explicitly, we write the spec- 
tral function in terms of symmetric and antisymmetric 
parts: ^(k,e) = As(k,e) + Aa{'k,£), where As(k,a;) = 
As(k, -w), Aa(k,cj) = -Aa{k, -uj) and uj = e-/i. Using 
these definitions, the c-axis thermopower may be written 



IV. APPLICATION TO Bi2212 

Each momentum state in the energy-dependent con- 
ductivity is weighted by the matrix element t^, and 
it is important to understand the effect that this has. 
In cuprates with a simple tetragonal structure and one 
Cu02 plane per unit cell, it is well known that the tun- 
neling matrix element depends strongly on the in-plane 
momentum, and vanishes along the zone diagonals 2^. 
The origin of this effect is that the hopping process occurs 
via the Cu 4s orbital. The overlap of this orbital with the 
in-plane Cu-0 hybrid orbitals has dj,2_j,2 symmetry. Be- 
cause the c-axis tunneling matrix element vanishes along 
certain directions in fc-space, one of its most important 
roles in c-axis transport is to select the region of /c-space 
which makes the dominant contribution [2^ , [25| . 

For Bi2212, two additional complications arise. The 
first is that the underlying lattice in this material is 
body-centered tetragonal. The off-set this creates be- 
tween neighboring Cu02 planes means that the hopping 
matrix element now vanishes along the zone-boundaries 
as well as along the zone-diagonals ^2^. 

The second complication is that Bi2212 is a bilayer 
compound, with two Cu02 planes per unit cell. A mini- 
mal description necessarily involves two hopping param- 
eters. One describes the hopping within the bilayer (and 
is responsible for the bilayer splitting), and the other in- 
volves hopping between different bilayers. In the limit 
where the bilayer hopping ti,i is much larger than the 
inter-bilayer hopping tc, the bilayer splitting can be sep- 
arated from the intercell hopping, giving the following 
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expression for the tunneling matrix element [27| : 

ik = zt4:tcCOs{kxa/2) cos{kya/2)x 

((cos fcajfl — cos A:j^a)^/4 + ao), (12) 

where a is the in-plane lattice spacing. The 'vertical hop- 
ping' parameter ao has been introduced phenomenologi- 
cally in Ref. [131 to account for hopping processes which 
are not assisted by Cu 4s orbitals. It is needed to fit 
the LDA band-structure, where a finite bilayer splitting 
is seen at the zone-center. Note that thermopower does 
not depend directly on the magnitude of tc, since it can- 
cels from the numerator and denominator of Eq. ([9]) . 

The spectral function in Bi2212 has been extremely 
well characterized by ARPES 2^. However, it is im- 
portant to remember that the ARPES intensity /(k, us) 
is proportional to /(aj)A(k, w), and so only the occu- 
pied states are probed. Norman et al (29l . [soj have pro- 
posed a technique to eliminate the effects of f{Lo) from 
ARPES data. Under the assumption that the spectral 
function is particle-hole symmetric at low energies and 
momenta close to the Fermi wavevector, the symmetrized 
intensity I{lo) + is just the spectral function (con- 

volved with the resolution). However, since the c-axis 
thermopower is non-zero (c/ Fig. [1]), we know that this 
assumption of particle-hole symmetry cannot be strictly 
correct. We can use the c-axis thermopower to obtain a 
measure of this asymmetry and to indicate whether im- 
portant physics is being lost by an analysis that assumes 
the spectral function to be symmetric. 

Let us consider the high temperature regime, i.e. at 
temperatures well above the pseudogap onset tempera- 
ture T* [sH. The temperature T* decreases with doping, 
and is eventually cut-off by Tc in the overdoped regime. 
Above T* , ARPES suggests that spectral function for the 
occupied states becomes only weakly dependent on tem- 
perature [2^, and hence the temperature dependence of 
the thermopower is dominated by the Fermi functions. 
Assuming that the energy-dependent conductivity varies 
slowly on the scale set by fc^T, the Sommerfeld approx- 
imation may be used to evaluate the thermopower: 

klT a'^jep) _ n^kskBT ^^^^ 

3 e cTcieF) 3 e eo 

The thermopower in this regime is therefore expected 
to be linear in temperature. In Fig. [1] we show the re- 
sults of fitting this form to the c-axis thermopower in the 
appropriate temperature range. The sublinear behavior 
at very high temperatures has two possible causes. One 
possibility is that the spectral function has some stronger 
temperature dependence not seen at lower temperatures. 
A second possibility is that at higher temperatures the 
Fermi functions explore crc(e) over a greater energy range. 
This may invalidate the Sommerfeld approximation, al- 
though the exact temperature dependence which results 
will depend sensitively on the exact shape of crc(e) [s^ . 

The quantity adep) / u'ciep) defines an energy scale cq- 
For free electrons in d-dimensions where the energy de- 
pendent conductivity varies as a power-law in energy (see 



Sec. |V]), this energy scale is proportional to the Fermi 
energy: eg = "itp/d. The extracted values of the energy 
scale Co are shown in the inset to Fig. [TJ The energy 
scale increases with the number of carriers p, while the 
sign reflects the sign of the charge carriers (holes) . 

We can obtain an independent measure of the expected 
scale of CT^/cTc from tunneling experiments on this mate- 
rial [33]. The tunneling conductance from an STM tip is 
related to the local density of states, — 2 ^(k, e). 
At high temperatures, the spectrum is fairly smooth and 
temperature-independent, but is not constant in energy. 
If we use the data of Ref. [3^ to estimate the value of the 
quantity v/v' over the band, we find a value of — 0.68eV 
for an underdoped (Tc = 83K) sample, and — 0.84eV for 
an overdoped (Tc = 74. 3K) sample. These energy scales 
are of the same order of magnitude, and have the same 
sign as those extracted from the c-axis thermopower. 

Further support for the assertion that the spectral 
function becomes only weakly temperature-dependent 
above T* comes from the c-axis resistivity. Using the 
Sommerfeld approximation, we would expect that in this 
temperature regime the c-axis resistivity should be ap- 
proximately independent of tem per ature. This is in good 
agreement with experiment [s^ l35j. where the c-axis re- 
sistivity becomes only weakly temperature dependent at 
high temperatures. 

We now turn to the behavior of the thermopower at 
lower temperatures. Here, the thermopower is seen to de- 
viate from linearity, pass through a maximum, and then 
finally fall to zero below Tc. Two effects become impor- 
tant in this regime. The first is the effect of the pseudo- 
gap, particularly in the underdoped and optimally doped 
samples. The effect of the pseudogap on the energy de- 
pendent conductivity is modulated by the tunneling ma- 
trix element, which selects the regions of fc-space which 
make the dominant contribution. Because vanishes 
along the zone-boundaries in Bi2212, the region around 
(tt, 0) point where the pseudogap is largest does not con- 
tribute significantly. Nevertheless, we would still expect 
to see some suppression of the energy-dependent conduc- 
tivity around the Fermi energy, although possibly to a 
lesser extent than that observed in the density of states. 

This suppression of the energy-dependent conductiv- 
ity has been linked to the increase in the c-axis resistiv- 
ity below T* [H, [131 • Since the conductivity appears as 
the denominator in the expression for the thermopower, 
it is tempting to associate the decrease in the conduc- 
tivity with the increase in the thermopower. The tem- 
perature at which the upturn in the thermopower occurs 
decreases with doping, which is consistent with the be- 
havior of the pseudogap opening temperature T*. This 
argument neglects the possible temperature dependence 
in the numerator of ([H]), which will more strongly reflect 
the changing asymmetry of the spectral function in the 
pseudogap regime. It would be interesting to combine 
both c-axis thermopower and resistivity data together to 
form the quantity acSc. This would indicate the extent 
to which the upturn in the thermopower is due to changes 
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FIG. 3: In-plane and out-of-plane thermopower of single 
crystals of the parent compound Bi2(Sr,La)2CaCu208 (after 
Ref. 

in the asymmetry of the spectral function. However ac 
has not been pubHshed for this crystal. 

Very close to Tc, a second possible contribution to the 
upturn in the thermopower are superconducting fluctu- 
ations. To our knowledge, the fluctuation contribution 
to the c-axis thermopower has not been calculated pre- 
viously, loffe et al. [3g| have calculated the fluctuation 
contribution to the c-axis conductivity. Near to Tc, they 
have identified a contribution which leads to a decrease 
in the conductivity, which is subsequently suppressed by 
the usual divergence at Tc- This decrease in the conduc- 
tivity suggests a possible explanation for the sharp rise 
in the c-axis thermopower close to Tc- This argument is 
particularly relevant for the p = 0.20 sample. The pseu- 
dogap, as determined from specific heat measurements 
[39|, is expected to vanish around p — 0.19, and hence 
may not explain the upturn in the thermopower in the 
p = 0.20 sample. 

For completeness, we conclude our discussion of Bi2212 
by considering the c-axis thermopower of the undoped 
parent compound. Fig. [3] shows Sc and Sab for the par- 
ent (r = 0) insulator Bi2(Sr,La)2CaCu208, reproduced 
from Ref. [17| . Remarkably, the thermopower of the par- 
ent compound is nearly isotropic, which suggests that 
the mechanism of transport in the insulating state is the 
same in- and out-of-plane. It has been proposed that 
the in-plane resistivity and thermopower [40[ can be ex- 
plained in terms of three-dimensional variable-range hop- 
ping [4l|. This approach gives Sab ^ Vt at low tem- 
peratures, and hence one may infer from the data that 
Sc VT also. Under the assumption that our theory 
of c-axis thermopower still applies to the parent com- 
pound, this result would provide information on the tem- 
perature dependence of the spectral function. However, 
this assumption may be invalidated by the fact that the 
parent compound is magnetically ordered, and so corre- 
lations between layers may have a significant effect on 
c-axis transport. 



In summary, by comparing our analysis of the c-axis 
thermopower to published experimental data on Bi2212, 
we can draw the following conclusions. By assuming 
that the spectral function becomes independent of tem- 
perature above T*, we have shown that the c-axis ther- 
mopower is expected to be linear in temperature and the 
c-axis resistivity temperature independent, at least while 
the Sommerfeld approximation remains valid. Compari- 
son with the experimental data has confirmed these pre- 
dictions, and so we conclude that both the symmetric and 
antisymmetric parts of spectral function become weakly 
temperature dependent above T*. The antisymmetric 
part of the spectral function is characterized by energy 
scales consistent with those seen in tunneling. Just above 
Tc, the thermopower shows an upturn whose microscopic 
origin could be related to the pseudogap physics and/or 
superconducting fluctuations. A comparison with resis- 
tivity measurements on the same crystals would reveal if 
this were due to changes in the antisymmetric part of the 
spectral fmiction. 

V. RELATION OF THERMOPOWER TO 
ENTROPY 

In this section we will discuss the relationship between 
thermopower and entropy in general terms. Recall that 
the thermopower or Seebeck coefficient of a sample is de- 
fined as the ratio of the induced voltage gradient to an ap- 
plied temperature gradient under conditions of zero elec- 
trical current flow, S — {dV/dT),j=Q. Under isothermal 
conditions the heat current density, jq, associated with 
a given electrical current density, j , is given by j'q = Hj , 
where H is the Peltier coefficient. The Kelvin/Onsager 
relations [l^. El show that the Seebeck and Peltier co- 
efficients are related by H = ST . It follows that we can 
write the thermopower as 

S = n/T={js/jhT=o, (14) 

where jg = jq/T is the entropy current density. The 
question is then how to interpret final equality on the 
right-hand side of Eq. (|14p ? The conventional interpre- 
tation is to say that "the thermopower is the transported 
entropy per charge carrier". However, given that a cur- 
rent density is a rate of fiow, an alternative interpretation 
of the thermopower is "the rate of change of transported 
entropy with number of charge carriers ". We also need to 
define the terms 'transported entropy' and 'charge car- 
rier': is the transported entropy the same as the equilib- 
rium entropy, and what is the relevant number of charge 
carriers? 

To further explore these issues, let us examine the case 
of non-interacting electrons in a single band, with an arbi- 
trary dispersion relation. The thermopower is calculated 
from the Boltzmann equation in the relaxation time ap- 
proximation Q. If the system is anisotropic, then the 
thermopower will generally be a 2nd rank tensor quan- 
tity: here, we choose to focus on just one of the diagonal 
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components, Sx- In addition, we consider the quantities 
-S/N^e, S/Nhe and -{1/ e){dS / dN^)Ty , where is 
the number of electrons and A^^i is the number of holes. 
Straightforward calculations provide: 



— = —kn 



= — kfi 



deu{e)[f\nf + {l-f)\n{l-f)\ 



dev{e)f 



dev[e)[f\nf + {l-f)\n{l-f)\ 



dev{e){l-f) 



(15) 



(16) 



The positive nature of S/N^ and S/Nh follows since 
r'(e) = iy{e) > 0. 

Let us now consider the two limits described in the in- 
troduction, where exact relationships between the ther- 
mopower and entropy hold. For non-interacting electrons 
with a quadratic dispersion in d-dimensions, it can be 
shown that ax{e) esc T(e)e'^/^ and j/(e) cx e''/^""'^, giving 
r(e) cx e*^/^. Hence if the relaxation time is taken to be 
an energy-independent constant, then ax{e) oc r(e), and 
Sx = —S/Nee. For non-interacting holes {i.e. with an 
inverted parabolic dispersion) , a similar analysis leads to 
Sx = S/Nue. 

The second limit is at temperatures much larger than 
the bandwidth. This limit is most easily understood by 
rewriting p8|) as 



eSx^F[ax{e% (^) = f^[.(e)], (17) 



where 



FW)] 



J deg{e)\ 




(- 


df\ 
de) 


J 


f deg{e)(^- 


de) 


1 



(18) 



and where axie) — 2e^ J^k ■'^^(k)'''(e(k))(5(e — e(k)) is the 
energy-dependent conductivity and z^(e) — 2^j^^(e — 
e(k)) is the density of states. These expressions reveal a 
fundamental difference between S/Ne and S/Nh on the 
one hand, and Sx and {dS / dNe)T.v on the other. Both 
S/Ne and S/Nh can be thought of as the sum of (posi- 
tive) contributions from electrons and holes. This means 
that S/Nf. and S/Nh are always positive quantities. In 
contrast, Sx and [dS / dNf.)T,v should be thought of as 
the difference of contributions from electrons and holes. 
These quantities can be either positive if holes dominate, 
negative if electrons dominate, or zero if the system is 
particle-hole symmetric. 

It is also possible to express S/N^ and S/Nh using 
the functional F. Integrating by parts in the numer- 
ator and denominator of (|15p and (flB)) . we find that 
S/Ne = F[T[e)] and S/Nh = F[N, - r(e)], where 
r(e) = jLoo '^{e')de' is the energy-dependent electron 
number, and N^ — r(oo) is the total number of states. 

Assuming that the function g{e) varies slowly on a scale 
set by fcsT around e = /j, ~ e^, the Sommerfeld approxi- 
mation may used to evaluate the integrals in p^ . giving 



F[g{e)] 



(19) 



where g' = dg/de. Hence F\g\ is expected to be lin- 
ear in temperature while the Sommerfeld approximation 
holds. The functional F\g] measures the relative asym- 
metry of the function g{t) about the Fermi energy. The 
thermopower and {dS / dNf/)T,v can either be positive or 
negative, depending on the signs of a'^iep) and v'{ep). 



FW) 



J deg{e)\ 






J 


f degie)[^ 


df\ 
de) 


1 



A* 
T' 



(20) 



For the functions v{e) and (Jx{e)^ the upper limit in the 
integral is the bandwidth, and so for sufficiently large T 
the first term tends to zero. Since the chemical potential 
is defined by ji/T = ~{dS/dNi,)E,Vi we recover Heike's 
formula: 



eSx 



dS 



idNe)T,V 



/ dS 

\dNe 



= fc R In 



E.V 



(21) 

where x — Ne/Ncew is the filling fraction. This result 
is not restricted to non-interacting systems, and can be 
shown to hold in general (for a system with no orbital 
deg rees of freedom) from an analysis of the Kubo formula 
|la |. Note that the same analysis does not apply to S/Ne, 
since the function r(e) is non-zero outside the bandwidth, 
and hence the first term in ([^D]) is finite in the large T 
limit. 

We now consider a particular band-structure away 
from the two limits discussed above. We take the 2d 
tight-binding model with an energy-independent relax- 
ation time as a simple, yet non-trivial example, for which 
e(k) = -2t(cos kx(i-\-cos kyd) . The relaxation time is also 
taken to be energy-independent: this crude approxima- 
tion is sufficient for our current purpose. 

Fig. Ufa) shows the density of states and energy- 
dependent conductivity for this model. In Fig. |ljb) 
we plot the thermopower, entropy per electron 
—S/NeC, entropy per hole S/Nhe and the quantity 
— {l/e){dS/dNe)T.v at a fixed temperature T = 0.2t. 
The thermopower is negative for x < 1 where electrons 
are the dominant carriers, and positive for a; > 1 where 
holes dominate. At half filling the system is particle-hole 
symmetric, and S = Q. At low filling where the dominant 
carriers are electrons, the thermopower tends to the en- 
tropy per electron —S /Nee. This is shown more clearly in 
Fig- Sic), and follows from the fact that the dispersion re- 
lation is parabolic close to k = (0, 0). At large filling, the 
dominant carriers are holes with an inverted parabolic 
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FIG. 4: Thermodynamic and transport quantities calculated 
for the 2d tight-binding model with an energy-independent 
relaxation time, (a) Density of states and energy-dependent 
conductivity, (b) Thermopower, entropy per electron, entropy 
per hole and derivative of entropy with respect to the number 
of electrons, plotted as a function of filling at fixed tempera- 
ture, (c) Ratios of the quantities presented in (b). 



dispersion, and the thermopower tends to the entropy 
per hole, S/NhC. Neither —S/N^e nor S/Nhe are good 
approximations to the thermopower close to half-filling. 

In contrast, —{l/e){dS/dNe)T.v has a qualitatively 
similar doping dependence to the thermopower, in the 
sense that it is negative for x < 1, positive for x > 1 
and passes through zero at half-filling where the sys- 
tem is particle-hole symmetric. The difference between 
S and —{1/ e){dS / dNe)T.y arises because each fc-state 



in the energy-dependent conductivity is weighted by the 
factor Wa;(k)^T(e(k)). This weighting factor is absent in 
the density of states. In essence it is this factor which 
converts entropy into 'transported entropy', and which 
converts the number of electrons into the number of 
'carriers'. Were it not for this factor, we would have 
S = —{l/e){dS/dNe)T.v for any band-structure, and 
for any temperature and filling. In the present exam- 
ple, this factor eliminates the strong features near the 
band edges and center which are present in the density 
of states, and which lead to the extra structure seen in 
-il/e){dS/dN,)T,v. 

The advantage of —{l/e){dS/dNe)T,v compared to 
the entropy per electron/hole emerges when we attempt 
to extend this analysis to interacting systems. Although 
the general expression for the entropy in an interacting 
system is extremely complicated [13, |43] , a relatively sim- 
ple expression may be derived for [dS / dNe)T.v ■ This 
quantity may be calculated by starting with the expres- 
sion for the total number of electrons: 



N,. 



2 /'de^A(k,e)/(e). (22) 



The derivative of this expression with respect to temper- 
ature at fixed particle number vanishes, since the total 
number of particles is conserved. Taking the derivative, 
and using the relation {dS/dNe)T.v — —{dfi/dT)]yy 
El, we find 



as 



T 



de) 



^^djjs) T ^ ( de 



(23) 

The first terms in the numerator and denominator come 
from the intrinsic temperature and chemical potential de- 
pendence of the spectral function. Note that the kinetic 
part of the spectral function does not contain the chem- 
ical potential, as this has been shifted into the Fermi 
function. Hence all temperature and chemical potential 
dependence comes from the self-energy. 

In summary, we have seen that even in a very 
simple example, the 2d tight-binding model with an 
energy-independent relaxation time, there is no uni- 
versal relationship between S and —S/Nee, S/Nhe or 
— {l/e){dS/dNe)T,v- Exact relations only exist in par- 
ticular hmits. We argued that [dS /dNf.)T,v may be 
the more appropriate quantity to compare with the ther- 
mopower, since both reflect the particle-hole asymmetry 
of the system. In a non-interacting system, the difference 
between these two quantities arises because {dS /dNe)T,v 
measures the relative asymmetry of the density of states, 
while the thermopower measures the relative asymmetry 
of the energy-dependent conductivity. A comparatively 
simple expression for [dS / dNe)T.v may also be obtained 
in an interacting system. 

In the next section, we will first make an experimental 
comparison between the c-axis thermopower of Bi2212 
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and the entropy per particle obtained from specific heat 
measurements, and then use Eq. (j23p to make a theoret- 
ical comparison with — {1 / e){dS / dNe)T.y ■ 



VI. RELATING THE C-AXIS THERMOPOWER 
TO THE ENTROPY 

Before we compare the c-axis thermopower with the 
entropy, we first offer an argument as to why one might 
expect a close relationship between the two. In the 
previous section, we compared the entropy for non- 
interacting electrons with the thermopower calculated 
from the Boltzmann equation in the relaxation time ap- 
proximation. We argued that one of the reasons why the 
thermopower and the entropy are different is that the 
relaxation time only enters the thermopower. This as- 
sertion is not quite correct. The scattering mechanism 
(impurities, phonons, electron-electron etc) gives rise to 
a self-energy, and will in general change the equilibrium 
properties of the system as well. One might imagine that 
if this self-energy were included in the calculation of the 
entropy, it would yield a closer relationship to the ther- 
mopower. In general this argument fails, since vertex cor- 
rections mean that the transport scattering rate is not the 
same as the imaginary part of the self-energy. However, 
this is not necessarily the case for c-axis transport, since 
we have argued in Sec. IIIII that vertex corrections may 
be neglected. Hence, c-axis thermopower provides an 
'optimal' case for making the comparison between ther- 
mopower and entropy. 

We now compare the experimentally measured ther- 
mopower with the equilibrium entropy obtained from 
specific heat measurements. For convenience, we re- 
produce in Fig. [Sfa) the thermopower data of Fig. [1] 
from to 200K. Fig. [DJb) shows the electron entropy 
per Cu (sy) reported in Ref. [HI, from which we calcu- 
late s^/eSc in Fig. [SJc). The monotonic increase in the 
entropy with temperature is very different from the in- 
plane thermopower [T7| . but is much closer to that of 
the c-axis. Now, if the thermopower measures the en- 
tropy per charge carrier (and the relevant charge carriers 
are the doped holes), then the ratio s^/eSc should be 
independent of temperature and equal to p, the number 
of charge carriers per Cu. Fig. [S^c) shows that the ra- 
tio is nearly independent of temperature above around 
lOOK, although deviates strongly as the transition is ap- 
proached. Also, although the ratio is clearly not equal to 
p, it does increase with increasing doping. The experi- 
mental results indicate that while the c-axis thermopower 
is certainly more closely related to the entropy than is the 
in-plane thermopower, there is at best only a qualitative 
correspondence between the two. 

We now turn to the theoretical comparison of the ther- 
mopower and —{l/e){dS/dNe)T,v- Eq- ((23)) bears a 
strong resemblance to the c-axis thermopower ([9]), al- 
though there are clearly some differences. The tunnel- 
ing matrix element appears in the thermopower but not 




Q\ , \ , 1 

100 200 

Temperature (K) 

FIG. 5: (a) c-axis thermopower Sc of Bi2212 single crystals, 
reproduced from Fig.[T] The data are plotted as eSc in units of 
ks, to allow direct comparison with the entropy, (b) electron 
entropy per Cu (s„) of Bi2212 single crystals measured by 
Loram et al. plotted in a unit of fes per Cu. (after Ref. ^H). 
(c) Sv/eSc calculated from (a) and (b). 



in {dS/dNe)T,v- This was to be expected, as we saw 
a similar difference in Sec. |V] where the group velocity 
entered the thermopower but not {dS / dNe)T,v ■ Since 
the tunneling matrix element depends strongly on the 
in-plane momentum, the dominant contributions to the 
thermopower and {dS / dNe)T.v niay come from different 
regions of fc-space. 

The first terms in the numerator and denominator of 
{dS / dNf,)T.v come from the intrinsic temperature and 
chemical potential dependence of the spectral function, 
and do not appear in the thermopower. At high temper- 
atures, we have argued that the spectral function may 
become rather temperature and doping independent, and 
hence these terms could be neglected. The same is not 
true at low temperatures where the spectral function is 
strongly temperature dependent and very sensitive to 
doping. 

The second terms in the numerator and denominator of 
{dS /dNe)T,v closely resemble their counterparts in the 
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thermopower, except that the spectral function appears 
hnearly rather than quadraticaUy. Note that if ^(k, e) oc 
A^{k,e), then because {dS/dNe)T,v is a quotient it will 
be unchanged under the replacement A(k, e) ^^(k, e). 
If the spectral function has any strong features, we may 
expect a significant difference between its first and second 
powers. For a spectral function which is rather broad and 
featureless, this difference will be less pronounced. 

For the particular case of Bi2212, the conclusion is 
that at high temperatures the spectral function is rather 
featureless, and the self-energy temperature and doping- 
independent. Hence, we may expect the thermopower 
and —{l/e){dS /dNe)T,v to have similar temperature de- 
pendence, although the coefficient of proportionality be- 
tween the two will not necessarily have any physical 
significance. At lower temperatures where the spectral 
function does depend strongly on temperature and dop- 
ing, we would not expect any kind of simple relation- 
ship between the thermopower and —{l/e){dS/dNe)T.v- 
Finally, it is interesting to note that the relationship 
we have predicted between the c-axis thermopower and 
— {l/e){dS/dNe)T,v is essentially identical to that be- 
tween the thermopower and the entropy per hole seen in 
Fig.Ii 

VII. SUMMARY 

In this paper we have investigated the c-axis ther- 
mopower in a quasi-two dimensional material. We have 
shown that within the tunneling Hamiltonian formalism 
the c-axis thermopower may always be written as a Mott 
formula in terms of the energy-dependent conductivity, 
crc(e) = 2e^ X)k *k^^(^' This implies that the c-axis 
thermopower is a measure of the relative asymmetry of 
the spectral function around the Fermi energy. 

Using this theory, we have shown that the c-axis ther- 
mopower and resistivity imply that the symmetric and 
antisymmetric parts of the spectral function become ap- 
proximately temperature independent at high tempera- 
tures. At lower temperatures, we have suggested that the 
up-turn in the c-axis thermopower may reflect the open- 
ing of the pseudogap. It is also possible that very close 



to Tc superconducting fluctuations may also contribute 
to this upturn. The main conclusion of our analysis is 
that the c-axis thermopower is an effective probe of the 
particle-hole asymmetry in the system. We believe that a 
systematic analysis of c-axis resistivity and thermopower 
data could provide valuable insights about the spectral 
function, particularly in the pseudogap regime. The c- 
axis resistivity and thermopower data together would 
also form a stringent test of any spectral function ob- 
tained phenomenologically, or from a microscopic model. 

The second theme of the paper is the relationship 
between thermopower and entropy. We discussed this 
relationship in general terms, and suggested that the 
derivative of the entropy with respect to particle number 
may be more suitable to compare with the thermopower, 
rather than the entropy per particle. We argued that the 
c-axis thermopower of a quasi-two dimensional system is 
a case where there might exist a close relationship with 
the equilibrium entropy. In the case of Bi2212, we showed 
experimentally that, at best, there is only a qualitative 
correspondence between the c-axis thermopower and the 
entropy per hole at moderate temperatures. Having de- 
rived expressions for both the c-axis thermopower and 
{dS / dNe)Ty , we were able to make an explicit com- 
parison between the two. Our analysis indicated that 
the relationship between these two quantities would be 
similar to that between the thermopower and entropy 
per hole seen experimentally. Because of the existence 
of a relatively simple expression for {dS / dNe)T.v , even 
in an interacting system, we believe that this quantity 
will prove useful in investigating the relationship between 
thermopower and entropy in other systems. 
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